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Abstract —This paper addresses the problem of rate selec¬ 
tion for the cooperative hybrid automatic repeat request with 
chase combination (HARQ-CC) system, where time correlated 
Nakagami-m fading channels are considered. To deal with this 
problem, the closed-form cumulative distribution function (CDF) 
for the combine SNRs through maximal ratio combining (MRC) 
is first derived as a generalized Fox’s H function. By using this 
result, outage probability and delay-limited throughput (DLT) are 
derived in closed forms, which then enables the rate selection 
for maximum DLT. These analytical results are validated via 
Monte Carlo simulations. The impacts of time correlation and 
channel fading-order parameter m upon outage probability, 
DLT and the optimal rate are investigated thoroughly. It is 
found that the system can achieve more diversity gain from less 
correlated channels, and the outage probability of cooperative 
HARQ-CC system decreases with the increase of m, and etc. 
Furthermore, the optimal rate increases with the number of 
retransmissions, while it decreases with the increase of the 
channel time correlation. 

Index Terms —Hybrid automatic repeat request, chase com¬ 
bination, time correlation, Nakagami-m fading channels, delay- 
limited throughput, rate selection. 

I. Introduction 

Hybrid automatic repeat request (HARQ) is a reliable and 
effective transmission technique to overcome the channel im¬ 
pairments without the knowledge of channel state information 
(CSI). In HARQ scheme, the transmission information is first 
encoded via forward error-correcting encoder. The encoded 
packet will be transmitted in a number of HARQ rounds until 
the receiver successfully decodes the packet. On the basis of 
whether the previously failed packets are used for decoding or 
not, HARQ scheme can be classified into two types, i.e. the 
Type I and the Type II HARQ. For the Type I, the previously 
failed packets will be discarded. On the contrary, for the Type 
II, the previously failed packets will be stored for reutilization. 
By using different decoding techniques to deal with all the 
received packets, the Type II can be further divided into two 
schemes, i.e. HARQ with chase combining (HARQ-CC) and 
HARQ with incremental redundancy (HARQ-IR). For HARQ- 
CC, the same packet will be transmitted in each HARQ round, 
and the technique of maximal rate combination (MRC) is 
adopted to combine all the received packets. However, new 
parity bits will be sent in each HARQ round for HARQ- 
IR, and code combining method is used for decoding the 
received packets. It is obvious that the Type II is superior 
to Type I because of exploiting the previously failed packets. 
Meanwhile, HARQ-IR generally performs better than HARQ- 


CC in terms of data transmission rate because of utilizing 
complex code combining approach. However, the HARQ-CC 
scheme has a low complexity, and a low hardware requirement, 
especially the buffer size. Thus our focus is put onto HARQ- 
CC in this paper. 

Performance of HARQ-CC has been frequently discussed 
in the literature. Generally, the performance evaluation of 
HARQ-CC depends on the distribution of the combining 
signal-to-noise ratio (SNR) through MRC. In ID, the opti¬ 
mal power assignment is presented for quasi-static Rayleigh 
fading channels. Because of the channels are constant in 
each HARQ round, the combine SNR follows exponential 
distribution. In 0, optimal power allocation is studied for 
i.i.d. Rayleigh fading channels. To simplify the analysis, a 
truncated probability density function (PDF) of the combine 
SNR was derived. The truncated PDF become more accurate 
for a higher transmitted SNR. Both the HARQ-IR and the 
HARQ-CC schemes are studied in cooperative system over 
i.i.d. Rayleigh fading channels in 0, where the exact PDF of 
the combine SNR is derived for evaluating HARQ-CC. 

In the practical scenarios, users with low-to-medium mobil¬ 
ity typically suffer from time-correlated fading El- An optimal 
rate adaption scheme is proposed for time-correlated Rayleigh 
fading channels, where the optimal rate is obtained through 
maximizing the delay-limited throughput (DLT) by using the 
Gaussian approximation for the combine SNR. In Q, another 
optimal rate selection scheme is presented for time-correlated 
Nakagami-m fading channels via minimizing the allocated 
resource. In that paper, the combine SNR is approximated as 
a Gamma random variable. An adaptive power allocation is 
proposed for time-correlated Rayleigh fading channels in El 
by using the same approximation for the combine SNR as El- 

Literatures seldom put their interest upon the cooperative 
HARQ systems, especially for the time-correlated channels. 
Due to the involvement of cooperative communications and the 
presence of time-correlation, the rate selection becomes more 
challenging than the prior works. In this paper, we investigate 
the rate selection problem in the cooperative HARQ-CC 
system over time-correlated Nakagami-m fading channels. The 
PDF of the combine SNR is first derived in closed-form as a 
generalized Fox’s H function. Based on the analytical result, 
outage probability and DLT are then derived, which then 
enables the rate selection by maximizing the DLT. The impacts 
of time-correlation, the maximum number of transmissions and 
channel fading-order parameter m on outage probability, DLT 
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Fig. 1. System model 


and the optimal rate are finally investigated for cooperative 
HARQ-CC system over time-correlated Nakagami-m fading 
channels. 

The remainder of the paper is organized as follows. In 
Section [II] the system model is presented. In Section [Till the 
distribution of the instantaneous combine SNR via MRC can 
be derived in terms of the generalized Fox’s H function. 
By applying these analytical expressions into performance 
analysis, the outage probability and DLT for HARQ-CC are 
derived in Section |IV] and an optimal rate selection scheme is 
presented. The Section Ivl draws some important conclusions. 


II. System Model 


In this paper, we consider a cooperative system of three 
nodes, as shown in Fig. Q] Prior to the transmission of b 
information bits to the destination, these bits are encoded and 
transmitted with a rate of 7 Z. Denote the transmitted signal 
as x. The power of transmitted signal x is assumed to be 
P s = K(x 2 ). The same signal will be transmitted in each 
HARQ round. The transmission procedure consists of two 
phase by employing HARQ-CC into cooperative system. For 
the phase I, the source broadcasts the transmitted signal x to 
the relay and the destination. The received signal corrupted by 
the block Nakagami-m fading channels at the destination and 
the relay in the fcth HARQ round can be expressed as 


Usd — hg D x + n SD ( 1 ) 

Vsr = hg R x + n k SR (2) 

where h k SD and h k SR represent the channel coefficients as¬ 
sociated with the links of source-to-destination and source- 
to-relay, respectively; n,g D and n RR are defined as additive 
white gaussian noises with mean zero and variance a 2 , i.e. 
n SD ( n %R ) ~ iV( 0 , er 2 ). Unlike the previous research upon 
HARQ-CC scheme Q, G), 0 , time-correlation of chan¬ 
nels is considered in this paper. The time-correlated chan¬ 
nels are modeled as a multivariate Nakagami-m distribution 
with exponential correlation. Specifically, the amplitude of 
channel gain \h k SD \ follows a Nakagami-m distribution, i.e., 
\hg D \ ~ Nakagami (m,rig D ). Hereby, the squared ampli¬ 
tude of channel gain \h k SD \ 2 are Gamma-distributed with pa¬ 
rameters (rn,£l k SD /m), i.e. \h k SD \ 2 ~ Q(m,Xl k SD /m). Here a 
widely used exponential time correlation model 0 is adopted 
and the cross-correlation between \hg D \ 2 and \h J SD \ 2 can be 
derived as 


PSD ~ 


cov^sp ] 2 ,\h J SD \ 2 ) 
/ \Jvar(\h i SD \ 2 )var(\h J SD \ 2 ) 


~ Kso tfsD ( 3 ) 


The same time-correlation model applies to \h RR \, that 
is, \hg R \ ~ Nakagami (rn,V,g R ), and the cross-correlation 

between \h l SR \ 2 and \h J SR \ 2 is denoted by p l RR = X l SR X J SR . 

In order to decode the received signals at the destination 
in this phase, MRC is adopted to combine all the previously 
received signals. Assume that M HARQ rounds are utilized, 
the combining signal is therefore formulated as 

tsd = w T y sd = w T h S dx + w T n S D ( 4 ) 

where y SD = {ysDfVsDr ■ ■ iV^dY defines the vector 
of received signals; w = {w\,W2, ■ • ■ ,WmY defines the 
weights of MRC; h SD = {hg D ,h 2 SD , ■■■ , h^ D } T and n SD = 
{rig D , rig D , ■ ■ ■ ,n^ D Y are defined as column vectors of 
channel coefficients and Gaussian noises in all HARQ rounds, 
respectively. Since the signal x has average power P s , the 
instantaneous output SNR for MRC after M HARQ rounds is 

_ M P s |w T h SD | 2 __ P s |w T h SD | 2 ^ 

Td, 1 r 9 'i 11 2 CH 

E||w T nsi)| 2 j (j2 ll w ll 

Since constants do not matter, one could always scale 
w such that ||w|| 2 = 1 . For MRC technique, the output 
SNR 7 m is maximized. By using Cauchy-Schwarz inequality, 
the expression 0 has a maximum only if w is linearly 

j^k * j. 

proportional to hsu, i.e. Wk = p^uf. Therefore, 7^ is 
rewritten as 



Similarly, the combine SNR at the relay after M HARQ 
rounds is formulated as 



Once the destination fails to decode the received signal y k R , 
the destination will feed back a negative acknowledgement 
(NACK) to the source node. The next HARQ round will be 
immediately triggered and the source retransmits the same 
signal x once receiving the NACK message. The number 
of HARQ-CC rounds is limited to K to avoid network 
congestion. Only if the relay successfully decodes the received 
signal, the transmission procedure will proceed to the phase 
II, that is, the source and the relay are cooperative to transmit 
the signal. In this phase, the received signal associated with 
the link between the relay and the destination is 

VrD = hRD X + n RD ( 8 ) 

where h k RD represents the channel coefficient associated with 
the link of relay-to-destination; n k RD is defined as an additive 
white Gaussian noise with mean zero and variance a 2 , i.e. 
n RD ~ N(0,cr 2 ). Similarly, \h RD \ ~ Nakagami (' m , Qrd)’ 
and the cross-correlation between \h RD \ 2 and \h RD \ 2 is de- 
noted by p l RD = X RD . Similar to the phase I, suppose 
that the relay successfully decodes the received signal after r 








HARQ rounds, then the resultant combine SNR relied upon 

{ ysD,y l RDi k e [1 ,M],l G [r + l,M]} is expressed as 
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Ps 
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SJt 


(9) 


It should be mentioned that h k SD , h k SR and h k RD are inde¬ 
pendent of each other. In the case of successively decoding the 
signal x at the destination no matter in the phase I or the phase 
II, an acknowledgement (ACK) message will be transmitted 
to the source. Consequently, the retransmissions for x will be 
ended, and the subsequent b information bits will be encoded 
for transmission in the same way. To address the problem of 
rate selection, the distribution of the combine SNR should be 
derived first. 


III. Closed-Form PDF and CDF for Instantaneous 
Combining SNRs 

From ©, 0 and 0, all the combine SNRs are expressed 
as the sum of correlated Gamma random variables (RVs). The 
nature of our work becomes to determine the distribution of 
sum of correlated Gamma random variables. Some previous 
literatures have studied this problem. In ( 8 ), 0, the PDF for 
the sum of correlated Gamma RVs is expressed as multi¬ 
fold infinite series. In 0, the PDF for sum of correlated 
exponential (a special case of Gamma) RVs is expressed as a 
confluent Lauricella hypergeometric function of multivariates, 
which is multi-fold integration. However, both the results 
obtained in 0 and 0 are hard to evaluate precisely because 
of involving of multiple folds. Imran S. A. in COl give a 
new result for the distribution of sum of independent Gamma 
RVs. Its PDF is derived as a generalized Fox’s H function. 
It is noteworthy that this alternative result presented in Col 
can be expressed in the form of a single-fold integration, 
which is readily computable. The authors in ITOl provide an 
efficient MATHEMATICA implementation of the generalized 
Fox’s H function, which can be evaluated fast and accurately. 
Following a similar approach as in m, we here aim to derive 
the distribution of the sum of correlated Gamma RVs. 

We first derive the distribution for Tdi = ZlLi z k> 
where z k = |#|^s,d| 2 - It can be proved that Zk follows 
Gamma distribution with parameter (m, VI k SD /to), i.e. Zk ~ 
G(m,Q RD ' /to), where Vl RD ' = Cl RD P s /a 2 . Without loss 
of generality, Zk is distributed as Zk ~ G(jn,0k), where 
0k = so /to. After some manipulations, it can be obtained 
that the cross-correlation between Zi and Zj is equal to p l RD . 
To cope with this problem, the following lemma is given for 
the PDF of Y m = Z£Li z k- 


Lemma 1 (PDF of the Sum of Gamma Random Variables 
with Exponential Correlation). The pdf of x can be derived 
in terms of the generalized Fox’s H function SUl Def A.57] 


as 

(y) = ff 4~ m x 

Td,i k=1 

(l - <5i -1 , l,m) , • • • , (l - <5 m - \ l,m) y 
(-4 - \ l,m), • • • , l,m) 

( 10 ) 

where {S k }kLi are defined as the eigenvalues of the matrix 
A = DC, where D is the M x M diagonal matrix with 
the diagonal entries {0k}ff—i, and C is the M x M positive 
definite matrix given by 

I_ V Psd 

c _ VPsd 1 

/ M, 1 / M, 2 

L V Psd V Psd 

Proof: Please refer to 0 and COl . ■ 

Accordingly, by using the result from Lemma |T] the cu¬ 
mulative distribution function (CDF) of Ym is given by the 
following Theorem. 

Theorem 1 (CDF of the Sum of Gamma RVs with Exponen¬ 
tial Correlation). For arbitrary m, the CDF of 1 can be 
derived in term of the Fox’s H function as ([72]) as shown on 
the top of next page. 

Proof: See Appendix [0 ■ 

There are some special cases for the expression (fl2l) as 
follows: 

1) In the case of M = 0, F ^ (y) = 1 holds since 

Td, i 

m 1 r e~ sy 

(y) = ±~ 7 T-f - ds = l-u(-y) (13) 

Td,i 27T* Jq S 

where u (y) represents Heaviside step function. 

2) For M = 1, because = Z\ follows Gamma 

distribution, it follows that 

Cl (») = Pr (a < y) = pAyT (m ,£) (14) 

where T (k,6) represents the lower imcomplete gamma 
distribution. 

3) If m is an integer, the generalized Fox’s H function 
reduces to Meijer-G function lfl2l Def. 9.30], thus 
F^m ( y) is rewritten in terms of Meijer-G function as 

Specifically, the corresponding CDF for Rayleigh 
fading channels can be obtained with m = 1 . 

In a similar way, the closed-form CDF F 1 ']} (y) for combine 

Tr 

SNR r y R can be derived. The remained work is to deter¬ 
mine the distribution of y'p f ■ using moment-generating 
function (MGF), the CDF of Td 2 ' s gi yen in the following 
theorem. Without loss of generality, we define yi = 
then 7 D j 2 = Zk=i z k + Lz=r+i Vu where z k and yi are 
independent of each other, and {yi}f^ 1 are exponentially- 
correlated Gamma RVs. It follows that yi ~ G(m, 0i), where 
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^ A 
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(15) 


9i = Pl l RD P s / (toct 2 ), and the cross-correlation between yt 
and y 7 is equal to p l RD . 

Theorem 2 (CDF of For arbitrary m, the CDF of 

yf'f can be expressed in term of the Fox’s H function as ( 1771 ) 
on next page, where {otk}ff— r . i are defined as the eigenvalues 
of the matrix B = FE, where F is the (M — r) x (M — r) 
diagonal matrix with the diagonal entries {9i}ffL r+1 , and E 


— r) x (M 

— r ) positive definite matrix given 

by 

1 

/ r +l,r+2 

V Prd 

/ r+l,M ~ 

V Prd 



/ ~r+2,r+l 

V Prd 

1 

/ r+2,M 

V Prd 


(16) 

/ M,r+1 

V Prd 

/ M,r+2 

V Prd 

i 




Proof: See Appendix iBl ■ 

IV. Optimal Rate Selection 


the relay successfully decodes the packet at round r. Therefore, 
to address the rate selection, the outage probability P R lt {r) 
and P’jj t (k\r) should be derived first. 

A. Outage Probability 

For P™ t (r), it is written as 

Pr 1 (r) = Pr (I r R < TZ) (20) 

where I R represents the total accumulated mutual information 
for MRC at the relay till M HARQ rounds, which is evaluated 
as 

Ir = log 2 (1 + Yr) (21) 

By substituting (l2ll into (l20t . it yields 

Pr* (r) = Pr (Yr < - l) = F™ [ 2 n - l) ( 22 ) 

The last step holds by applying the Theorem Q] 

For Pff t (k\r), it is given by 


An optimal rate selection scheme is proposed for coopera¬ 
tive HARQ-CC system over time-correlated fading channels 
via maximizing the DLT. Given the maximum allowable 
number of transmissions for a single packet K and data 
transmission rate 1Z, the DLT is given by m 

Tk = J2j (PS^fc - 1) - P D*(k)) (18) 

fc =i 

where P ^(k) defines the outage probability that the destina¬ 
tion fails to decode the received signals after k HARQ round. 
Specifically, P™ 4 (0) = 1, since it is impossible to decode 
the packet at round 0. By using the law of total probability, 
Pff t (k) is formulated as 

OO 

pS“(fc) = E ( P ^( r - !) - P °R Ut (r)) P° D ut (k\r), 1 < k < K 

r =1 

(19) 

where P™*(r) defines the outage probability at the relay 
after r HARQ rounds, and P™ ,t (fc|r) represents the outage 
probability at the destination after k HARQ rounds given that 


P™\k\r) = Pr (lf r < Uj (23) 


where iff represents the total accumulated mutual informa¬ 
tion for MRC at the destination till M HARQ rounds given 
that the relay successfully decodes the packet at the rth round, 
and is formulated as 


7-fc ,r 

L D 


log 2 (1 +7 d,i) p > k 
log 2 (l + 7D, r 2 ) . else 


(24) 


By putting (l24l > into < 1 2 3 b and then applying the Theorems 
|T]and [2] it follows 


P° D ut (k\r) 


(2 K -l),r>k 
F (2 J r (2 n - 1), else 


Hereby, by substituting d23l > and (l25l ) into (fl9l> . P'jf lk) can 
be expressed as <[26t on next page. 

In order to study the impact of channel correlation and 
fading severity on cooperative HARQ-CC system, we consider 
a common scenario for numerical analysis. Assume yr = 


























M M 

F% (y) = 1 - II 4“ m II — x 

^ D > 2 fc=l /=r+l 



M—bracketed terms (M—r) — bracketed terms 

' 

frh M 

n M+l,M+l 

(l — <5i , 1,to) , • • • , (l - 5m \ l,m), (1 - a r+1 ' , l,m) , ■ • • , (l - a M \ 

(0,1) , (-<5i -1 , 1 , to ) , • • • , (-5m' 1 , 1 ,m), (-a r+ i -1 , l,m) , • • • , (~a M ~ 

1,to),(1,1) 

1 ,1, to) 

e y 


M—bracketed terms (M—r) — bracketed terms 



( 17 ) 


k -1 


D 


! (fc) = 


r=l 


^ (2^-l)-F«(2^-l)) 
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(26) 




Fig. 2. Outage probability P^^K) vs. transmit SNR 7y. 


P s /a 2 , which is termed as transmit SNR. The mean power 
of each channel coefficient is assumed to be E(|/t<j£)|) = 0.5, 
= ^(\^rd\) = 1- Under these assumptions, the 
outage probability is plotted versus transmit SNR 

7 t for different p, m and K by setting R = 2 bps/Hz in 
Fig. [2] The presented results show a perfect match between 
the analytical and simulation results for p = 0 . 2 , 0.8 and 
to = 1, 2. From this figure, it can be observed that the outage 
probability decreases as transmit SNR 777 maximum trans¬ 
missions K and fading-order parameter to increase. Moreover, 
the outage probability increases as the correlation coefficient 
p increases. 

B. Optimal Rate 

To find the optimal rate, the relation between DLT T and 
transmission rate TZ is investigated first. An example is given 
in Fig. [3] by setting m = 2 and 77 - = 0 dB. From this figure, 
it can be observed that each curve has only one maximum 
point. This sole maximum point can be effectively obtained 
by adopting Golden section search. Hereby, for transmitting 
certain b bits information, the DLT is maximized by selecting 
an optimal rate which is determined from the following 


Fig. 3. DLT T vs. transmission rate 7 Z. 

optimization problem: 

maximize 7m 
n 

subject to TZ > 0. 

As exhibited in Fig. Q] for a fixed value of to = 2, 
the optimal rate for different K and p is plotted against 
transmit SNR. It can shown that the optimal rate increases 
as 7 t increases. For instance, the optimal rate increases from 
1 bps/Hz to 2.2 bps/Hz when 7 t varies from 0 to 10 dB 
for p = 0.2 and K = 2. If the maximum transmissions K 
is increased, the optimal rate increases as well. Moreover, the 
optimal rate decreases with the increase of time-correlation of 
channels. It is due to the fact that the channels with high time- 
correlation will degrade the time diversity gain of HARQ. 

V. Conclusions 

In this paper, an optimal rate selection scheme has been pro¬ 
posed for cooperative HARQ-CC systems, where the general 
time-correlated Nakagami-m fading model is considered. The 
involvement of cooperative communication and the presence 
of time-correlation make the problem challenging. The nature 
of rate selection is to choose the appropriate transmission 
rate to maximize the DLT. To cope with this problem, it is 
necessary to determine the distribution of the combine SNRs 





































Appendix B 
Proof of Theorem[2] 

Due to the independence of yfk=i z k and Ei! r+ i yi ’ the 
MGF of 7 ^ is given by 


M M,( S ) = nf =1 S k - m (6 k -' +s)- m x 

UZr + i^r m (ar 1 + s)- m 

The PDF of "ff ’2 can be straightforward derived as a 
generalized Fox’s H function via inverse Laplace transform. 
Furthermore, the CDF of 7 ^f’ 2 T is derived as ( fTTb by using the 
same method from Theorem [Q 
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Fig. 4. Optimal transmission rate 7 Z vs. transmit SNR 7 j-. 


via MRC. By means of the generalized Fox’s H function, 
the PDF and CDF of combine SNRs have been derived, with 
which the outage probability and DLT are then obtained. 

These analytical results have been validated by Monte Carlo 
simulations. It has been found from the results that the system 
can achieve more diversity gain from the less correlated chan¬ 
nels, and the outage probability decreases with the increase 
of fading-order parameter m, and etc. It is noteworthy that 
the metrics for evaluating the performance of the cooperative 
HARQ-CC system are derived in terms of the generalized 
Fox’s H function, which has a low computational complexity 
and can be evaluated accurately. Finally, the optimal rate to 
maximize the DLT has been presented. 


Appendix A 
Proof of Theorem!]] 


By integrating the PDF expressed in (flOb from 0 through 
y, the CDF of r y^ 1 is obtained, more precisely. 
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—m 1 x _ fc=i _ 

_ _ M 

fe=i r(i-s) n r">(i+{ fc - 1 +8) 

fc = l 


Sy ds 


e sy ds 


(27) 

After identifying the Mellin-Barnes contour integral in the 
above equation with the generalized Fox’s H function, the 
CDF of 7^7 is obtained as (fT2l i. 
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